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PREFACE 

This  paper  presents  the  results  of  a  study  of  Constant 
Absolute  Vorticity  Trajectories  applied  to  isentropic  surfaces. 
The  objectives  of  this  study  were:  first,  to  show  that  Constant 
Absolute  Vorticity  Trajectories  are  theoretically  better  applied 
to  forecasting  future  positions  of  parcels  on  isentropic  surfaces 
than  on  constant  pressure  surfaces  as  they  are  now  applied;  sec- 
ondly, to  show  quantitatively  the  actual  iniproveraent  of  the  tra- 
jectories on  isentropic  surfaces  by  comparing  statistically 
results  of  actual  forecasts  from  both  types  of  charts;  thirdly, 
to  suggest  how  the  improved  technique  can  aid  forecasting,  par- 
ticularily  upper  air  prognosis. 

This  work  was  undertaken  as  the  thesis  requirement  for  the 
degree  of  >'aster  of  Science  in  Aerology  at  the  U.  S.  Kaval  Post- 
graduate School,  Xonterey,  California,  during  the  academic 
year  1952-1953. 

The  author  is  indebted  to  Associate  Professor  George  J. 
Haltiner  of  the  Aerology  Department  for  his  very  valuable  assis- 
tance and  constructive  criticism  during  the  investigation.  He 
also  wishes  to  acknowledge  the  assistance  of  Professor  A.  Boyd 
Mewborn  of  the  Mathematics  Department  and  the  author ^s  wife, 
Pauline,  v;ho  contributed  so  much  help  with  the  laborious  task 
of  decoding,  computing,  and  entering  isentropic  charts. 
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I .   INTRODUCTION 

The  meteorologist  is  often  called  upon  to  forecast  for  periods 
in  excess  of  tvventy-four  hotirs  for  which  simple  extrapolation  is 
generallj''  insufficient.  One  raethod  of  forecasting  entails  prepa- 
ration of  prognostic  surface  charts  and  preparation  of  forecasts 
from  these  charts.  As  most  meteorlogists  Icnow,  laaking  accurate 
prognostic  surface  charts  is  easier  said  than  done.  One  of  the 
important  tools  used  in  estimating  the  position  and  intensity  of 
surface  systems  is  a  prognostic  500  millibar  chart.  At  present, 
much  attention  is  being  paid  to  the  problem  of  upper  air  prognosis, 
and  considerable  information  has  been  published  concerning  fore- 
casting the  500  millibar  surface,  the  most  recent  being  Forecasting 
in  the  Middle  Latitudes  by  Riehl  and  collaborators.  It  is  hoped 
that  this  paper  may  shed  some  light  on  a  facet  of  this  problem  and 
may  contribute  to  better  prognosis  and  forecasting. 

One  aid  to  prognosis  of  upper  air  charts  is  the  construction 
of  forecast  air  parcel  trajectories  based  on  the  principle  of  con- 
servatism of  the  vertical  component  of  absolute  vorticity,^'P +-f  j  , 
where  ^f  is  the  vertical  component  of  relative  vorticity  and  -f 
is  the  vertical  component  of  the  coriolis  force  known  as  the 
coriolis  parameter.  Rossby  Ts  |  initiated  the  vorticity  concept 
by  shelving  that  under  certain  restrictive  assumptions,  the  follow- 
ing relationship  holds: 


C^+f)  =  co>istdht. 


(1) 


The  assumptions  referred  to  above  are  as  follows: 

a.  The  atinosphere  is  baro tropic. 

b.  The  atmosphere  is  a  homogeneous,  incompressible  fluid. 

c.  Motion  is  purely  horizontal. 

d.  Friction  forces  are  neglected, 

e.  There  is  no  horizontal  divergence. 

Rossby  19,  pp  268-289]  then  integrated  his  fundamental  equation 
and  was  able  to  obtain  trajectories  for  particles  vdiich  conserve 
their  vertical  component  of  absolute  vorticity.  This  trajectory 
is  Imown  as  the  Constant  Vertical  Component  of  Absolute  Vorticity 
Trajectory  and  will  be  referred  to  as  a  CAV  trajectory. 

The  assumptions  necessary  to  derive  the  vorticity  equation 
according  to  Rossby  are  obviously  quite  restrictive.  Starr,  as 
editor  of  the  Journal  of  Meteorology  (June,  1945),  suggested  that 
CAV  trajectories  might  better  be  depicted  on  an  isentropic  chart 
than  on  a  constant  level  chart  as  there  are  no  solenoids  on  an 
isentropic  chart.  Namias  16,  pp  372-374/  states  that  from  obser- 
vations available  isentropic  surfaces  appear,  to  the  first  approxi- 
mation, to  be  substantial  surfaces,  i.e.,  surfaces  that  contain  the 
same  air  particles  from  day  to  day.  Influences  of  non-adiabatic 
nature  are  appreciable  over  long  periods  of  time,  but  are  usually 
insufficient  to  disrupt  the  fundamental  isentropic  flow  from  one 
day  to  the  next.  The  derivation  of  the  vorticity  equation  for  floiv 
on  an  isentropic  surface  requires  fewer  assumptions.  Such  restrictions 
as  purely  horizontal  flow,  a  barotropic  atmosphere,  and  a  homogeneous 


(2) 


incompressible  atmosphere  may  be  removed.  Now,  an  isentropic 
surface  is  a  surface  of  constant  potential  temperature.  In 
meteorlogical  processes,  potential  temperature  is  one  of  the  more 
conservative  elements  by  which  parcels  of  air  may  be  identified 
from  time  to  time.  Since  the  atmosphere  is  normally  stable, 
potential  temperatvire  increases  with  height  and  the  atmosphere 
may  be  considered  to  consist  of  an  infinite  number  of  isentropic 
surfaces.  Sir  Napier  Shaw  jlO,  p  263]  first  suggested  that  charts 
of  isentropic  svurfaces  should  be  dra^vn  as  the  motion  of  air  is  best 
resolved  on  isentropic  surfaces.  It  will  be  shown  in  this  paper 
that  the  prognostic  trajectory  based  on  the  conservatism  of  the 
vertical  component  of  absolute  vorticity  is  theoretically  and 
practically  better  suited  to  use  on  isentropic  charts  than  on 
constant  pressure  charts  as  it  is  now  applied. 
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II.  THEORETICAL  IrTVTISTIGATION 

The  motion  of  particles  on  a  substantial  surface  can  be 
described  with  Lagrangian  methods.  The  form  of  the  vorticity 
equation  for  flow  on  an  isentropic  surface  derived  here  using  the 
Lagrajigian  method  of  attack  is  due  to  iraltirer'"'.  LetX^M^i  and-fc 
be  normal  Cartesian  coordinates  and  let^X"  Y  ^  ^^^  ~Y      ^® 
Lagrangian  coordinates  where  }C   and  X.  are  coordinates  of  the 
horizontal  projection  of  the  isentropic  siu'face,  "£    ~  tL       ard  "V 
refers  to  time  for  the  isentropic  surface.  Then 


f/>§i 


It  follows  from  the  above  equations  that 

where  f-  eT -t<l'^    '  ^^^  strear^i  function  for  an  isentropic  surface. 

P     9 

In  a  similar  way  we  find 

^i.  <?X    <£)2-S'       c?)C      «D2:    c)«  e>X  (3) 


-'^Associate  Professor  G.  J.  Haltiner,  U.  S.  Naval  Postgraduate  School 

(4) 


Further,  since  the  motion  is  assur.ied  to  be  isentropic, 


ur- =  ^Zi  -f  a^  -y--u-^g-^  (5) 

c)>     OX     aT"  . 


Substituting  in  the  equations  of  motion 


(6) 


^     5^    ^^     S^     Z'  ^J 


the  expressions  1,  2,  3,  4,  and  5,  we  obtain 


c)V    c>X     <^X      c)X 


^v^  -f-  u.  c>v-'  ^-Lf  ^vr  —  —  ^  F    ^  ^J^ 

^>^        <^x        £>T  ax' 

Differentiating  (7)  with  respect  to  T  and  (8)  with  respect  to^ 


(7) 


(8) 


and  subtracting,  we  obtain 


Collecting  terns  _. 

^yL-^©J    L^x    c)xJl^t    <Dxj       ^ix:pT    <^xj 

.  ^  D    pec  ^  ^x>,l    __   ex.  c>f    +-x^J£   ^-ff^   -i-^^l 

^[5t    c)xJ  -     c)s:       c)X     L^r    c>tJ 


(5) 


and  introducing  vector  notation  we  obtain 


-{M^e]+'^■'^^.f.yl\.^=^■'^^o^^iV,,V4^ 


^irhere  the  operator  '^    ~  ^  ^ ,    n  g) 

Finally,  we  may  write  this  as 


A. 

cKt 


l^+fj+^ra+^j^^-V-O  (^) 


since'V=^t  for  any  particle.  Now,  assume  the  divergence  ten..^  .V=0. 
Then  we  may  write  the  vorticity  equation  in  the  familiar  form 


K^e+^J-^        or(^f^,4y^--f-^ 


(10) 


VJTiile  this  form  is  similar  to  the  form  derived  by  Rossby,  there 
are  some  fundamental  differences.  Parcel  motion  is  not  limited 
to  the  horizontal,  but  is  three  dimensional.  Moreover,  a  barotropic 
and  incompressible  atmosphere  is  not  required. 
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III.    techin'ique  of  investigation 

The  relative  vorticity  may  be  written  in  a  form  similar  to 
the  usual  expression  for  the  vertical  component  of  vorticity  in 
natural  coordinates 


^   R3   ON 


xirhere  P   and  ^  are  measured  on  the  isentropic  map.  If  we  select 
a  parcel  in  a  broad  uniform  flow  or  in  the  axis  of  a  jet  stream,  iX—o. 
P'urther,  if  the  parcel  is  at  an  inflection  point  in  the  f low,  ^  sr  O 
Under  these  conditions  the  CAV  trajectory  can  be  computed.  It 
will  be  noted  that  R   defined  above  is  not  exactly  equal  to  the 
projected  radius  of  curvature  of  the  streax.iline  on  the  isentropic 
surface.  However,  the  approximation  involved  here  is  normally 
less  than  those  iTiade  in  the  actual  computation  of  CAV  trajectories 
as  described  by  Fultz  [4,  p  13 J  and  revised  by  I7obus  ri2J  .  h'hile 
approximations  have  to  be  made,  it  must  be  remembered  that  the  true 
worth  of  a  technique  lies  in  the  statistically  proven  results 
obtained , 

The  simplest  method  of  obtaining  a  CAV  trajectory  is  through 
the  use  of  a  mechanical  device  such  as  the  so-called  "Wiggle-Wagon" 
devised  by  V/obus  [[121  and  currently  in  use  in  the  WBAN  Analysis 
Center,  V/ashington,  D.  C.  Original  studies  of  the  CAV  trajectory 
assumed  a  plane  earth.  This  assumption  entailed  errors  due  to 
distortion  when  compared  to  trajectories  along  a  spherical  earth. 
The  Wobus  "Wiggle-Wagon"  corrects  this  defect.  To  aid  meteorolo- 
gists who  do  not  have  access  to  the  "\\riggle-Wagon",  the  U.  3.  Kavy 
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Biireau  of  Aeronautics  Project  AROWA  i_llj  has  published  a  table 
for  computing  CAV  trajectories  raade  from  trajectories  traced  by 
the  "¥iggle-V/agon" .  This  table  will  be  used  to  compute  GAV  tra- 
jectories in  this  paper.  For  application,  it  is  necessary  to  Icnow 
speed,  direction,  and  the  latitude  of  the  parcel,  all  at  an  inflection 
point  in  the  flow. 

The  next  step  in  the  investigation  vi^as  to  compute  trajectories 
on  isentropic  and  500  millibar  charts  and  test  the  results  statisti- 
cally. A  rather  complete  compilation  of  upper  air  and  surface  data 
are  available  in  the  Northern  Haiiisphere  Surface  and  500  Ilillibar 
Charts  Series  published  by  the  U.  S.  V/eather  Bureau.  These  remarkable 
publications  not  only  have  a  daily  series  of  logically  analysed  surface 
and  500  millibar  charts,  but  also  contain  a  compilation  of  all  data 
used  to  plot  these  charts.  Data  were  taken  from  the  period  of 
January-February  1949.  From  this  data,  a  set  of  thirty  isentropic 
charts  was  plotted.  The  303K  isentropic  surface  was  selected  because: 
(1)  this  surface  was  often  close  to  the  so-called  "level  of  non- 
divergence"  described  by  Bjerknes  Flj  ;  (2)  this  surface  often  con- 
tained much  of  the  500  millibar  injid  speed  maximum  axis  (jet); 
(3)  the  303i:  isentropic  surface  was  usually  sufficiently  far  from 
the  surface  to  remove  effects  of  surface  friction  and  turbiJ-ence. 
It  was  desired  to  select  an  isentropic  surface  containing  much  of 
the  500  millibar  isotach  m^axixiuri  axis  so  computation  of  CAV  tra- 
jectories could  aid  in  the  prognosis  of  the  500  Eiillibar  surface. 
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Inasmuch  as  the  derivation  of  GAV  trajectories  on  an  isentropic 
surface  assumes  the  divergence  tenn  to  be  zero,  the  selection  of 
an  isentropic  surface  close  to  the  level  of  no  horizontal  divergence 
may  be  desirable.  The  effect  of  the  divergence  term  is  discussed  in 
Chapter  IV. 

After  the  heights  of  the  303K  isentropic  surface  were  plotted, 
all  available  wind  data  pertinent  to  the  surface  were  plotted,  Miere 
\dnd   observations  v/ere  available  from  stations  other  than  radiosonde 
stations,  an  interpolated  value  of  the  height  of  the  isentropic  surface 
and  the  wind  nearest  that  height  were  entered.  Because  of  the  fairly- 
dense  network  of  upper  air  \>dnd  observing  stations  in  the  United 
States,  sufficiently  accurate  streamline  analysis  was  possible  with- 
out the  time-consuming  computations  of  stream  function  values. 

Streamline  analyses  were  made  for  0300Z  each  day  in  the  period 
from  6  January,  1949,  to  22  January,  1949,  and  from  10  February,  1949, 
to  20  February,  1949,  For  each  chart,  a  number  of  CAV  trajectories 
were  computed.  The  following  conditions  had  to  be  met  before  choosing 
the  initial  point:  (1)  the  point  had  to  be  in  or  very  near  an  inflection 
point  in  the  flow  (undergoing  no  effects  of  curvature);  (2)  the  point 
had  to  be  in  or  near  the  axis  of  the  isotach  maximum  (jet  streani); 
(3)  the  point  must  be  at  an  actual  point  of  \id.nd  observation  to  obtain 
correct  initial  wind  speed  and  direction.  This  procedure  netted  from 
one  to  foiu:  points  suitable  for  computation  of  CAV  trajectories  per 
chart. 
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For  comparison  with  each  CAV  trajectory  plotted  on  the 
isentropic  chart,  a  companion  CAV  trajectory  was  computed  on  the 
500  niillibar  chart  of  the  sane  time.  It  was  not  too  difficult 
to  find  a  companion  inflection  point  on  the  500  millibar  chart  as 
the  isentropic  surface  was  usually  close  to  the  500  millibar  surface. 
After  the  CAV  trajectories  were  constructed  on  each  chart,  forecast 
wind  velocities  were  determined  by  moving  the  parcel  along  the  tra- 
jectory a  distance  equal  to  the  initial  ^mid  speed  multiplied  by 
the  number  of  hours  to  forecast  time.  The  forecast  times  selected 
were  twenty-four,  forty-eight,  and  seventy-t^iro  hours.  These  fore- 
cast ^mids  were  then  plotted  on  the  verifying  charts  and  verified. 
On  the  isentropic  chart,  verification  depended  on  the  forecast 
position  falling  within  the  observational  network.  If  the  forecast 
position  fell  within  the  network,  but  not  on  a  particular  station, 
the  forecast  v/inds  were  verified  by  linear  interpolation  between 
nearby  stations.  However,  care  was  taken  to  use  only  stations 
where  linear  interpolation  would  be  reasonable.  On  the  500  r.iillibar 
chart,  forecast  \^dLnds  were  verified  in  a  sir:iilar  manner  if  the  fore- 
cast position  fell  withj.n  the  observational  network  or  on  a  reporting 
station.  If  the  forecast  position  on  the  500  millibar  chart  fell  on 
neither  of  these  places,  then  contours  were  assumed  to  be  streamlines, 
and  gradient  x^rinds  were  measured  to  verify.  Since  verification  on 
isentropic  charts  depended  on  forecast  positions  falling  in  the  North 
American  network  or  on  a  particular  station,  the  number  of  possible 
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verifications  fell  off  drastically  after  forty-eight  hours.  However, 
vd.th  hemispherical  500  millibar  charts  available,  forecasts  made  on 
the  500  millibar  chart  could  always  be  verified.  All  500  millibar 
CAV  trajectories  were  carried  out  to  seventy-two  hours. 

In  verif3d.ng  forecast  \id.nds,  observed  wind  direction  error 
was  arbitrarily  selected  negative  if  the  observed  wind  direction  was 
to  the  right  (cloclcwise)  of  the  forecast  \rind   direction  and  positive 
if  the  observed  wind  direction  was  to  the  left  (coiuiter-clockwise) 
of  the  forecast  wind  direction.  Similarily,  wind  speed  error  was 
selected  negative  if  observed  wind  speed  was  less  than  forecast  and 
positive  if  the  wind  speed  was  greater  than  forecast.  All  forecast 
winds  and  verifications  were  tabulated  in  a  form  suitable  for 
statistical  testing. 

Kistograras  of  wind  speed  error  and  wind  direction  error  were 
prepared  for  twenty-four  and  forty-eight  hour  forecast  times  for  both 
the  isentropic  and  500  millibar  charts.  These  histograms  are  depicted 
as  Figures  1  through  4.  Wind  speed  errors  were  grouped  into  cells 
each  having  a  cell  interval  of  10  Imots  error.  The  cell  midpoints 
are  foriO, i'9;"^20  etc.,  laiots  error.  Since  the  cell  boundaries 
are  exactly  halfway  between  the  cell  midpoints,  some  of  the  obser- 
vations fell  on  the  cell  boundaries.  In  this  case,  it  was  deemed 
best  to  divide  such  observations  between  each  adjacent  cell.  Less 
trouble  was  encountered  in  grouping  wind  direction  error.  All  wind 
direction  errors  were  in  increments  of  10  degrees  and  values  of  0,±  10 
^ZO ,   etc.  degrees  of  error  were  taken  as  cell  midpoints. 
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After  the  data  were  grouped  in  histograms,  normal  curves  were 
fitted  in  accordance  with  the  method  described  by  Hoel  ([s,  pp  191-194J. 
After  obtaining  these  fitted  normal  curves,^  tests  were  performed  to 
determine  if  each  sajnple  of  data  could  reasonably  be  from  a  normal 
population.  In  all  cases  the  sample  distributions  could  reasonably 
be  said  to  be  from  norriial  populations,  although  in  two  cases  the 
value  of^   was  close  to  the  critical  value.  Restrictions  on  theV 
test  were  met  by  grouping  cells  with  too  few  frequencies  together. 
After  showing  normality  of  the  data,  the  distributions  of  data  from 
the  isentropic  charts  were  compared  to  similar  distributions  from 
the  500  millibar  charts  by  the  use  of  the  "F"  test.  The  "F"  test 
was  made  to  determine  if  compared  data  distributions  from  the 
isentropic  and  500  millibar  charts  could  statistically  be  from  the 
same  parent  normal  population.  The  purpose  here  is  to  show  that  a 
given  distribution  of  data  from  the  isentropic  charts  is  not  from  the 
same  normal  population  as  the  similar  data  distribution  from  the  500 
millibar  charts.  If  there  is  a  significant  difference,  then  inspection 
and  intuition  will  reveal  which  of  the  two  types  of  charts  are  better 
siiited  to  forecasting  by  means  of  CAV  trajectories. 
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IV.  RESULTS 

It  is  well  to  look  first  at  the  data  in  "everyday"  terms  to 
find  what  kind  of  resiiLts  a  forecaster  inay  expect.  If  it  is  arbi- 
trarily decided  that  a  wind  direction  forecast  will  verify  as  a  "hit" 
if  it  is  within   20  degrees  of  observed  direction,  then  we  find  that 
96.4^  of  the  twenty-four  hour  wind  direction  forecasts  from  CAV  tra- 
jectories on  isentropic  charts  verified  as  "hits".  By  comparison, 
76. 3, J  of  sirdlar  forecasts  from  500  millibar  charts  verified  as  "hits". 
With  forty-eight  hour  wind  direction  forecasts,  81,^  from  isentropic 
charts  verified  as  "hits"  and  41/?  from  500  millibar  charts. 

If  it  is  arbitrarily  decided  that  15  knots  of  wind  speed  error 
is  not  too  much  to  verify  as  a  "hit",  we  find  that  76,o  of  twenty-four 
hour  wind  speed  forecasts  from  trajectories  on  isentropic  charts 
verified  as  "hits".  Of  the  data  from  the  500  millibar  charts,  42.4/a 
of  twenty-four  hour  wind  speed  forecasts  verified  as  "hits".  Similarly, 
forty-eight  hour  vrlnd.   speed  forecasts  verified  78.5/b  from  isentropic 
chart  data  and  39%   from  500  millibar  chart  data.  So,  percentagewise, 
forecasting  wind  vectors  by  CAV  trajectories  on  isentropic  charts 
was  superior  to  forecasting  wind  vectors  by  GAV  trajectories  on  500 
millibar  charts  during  the  period  of  this  study. 

Figures  1  tlirough  4  are  histograms  of  error  betw^een  observed 
and  forecast  wind  speed  and  direction.  Applying  standard  statistical 
methods,  normal  curves  were  fitted  to  the  data.  The  ^  test  was 
applied  to  determine  if  the  distributions  of  error  were  likely  to  be 
from  a  normal  population.  Only  the  twenty-four  and  forty-eight  hour 
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errors  were  compared  and  tested  as  the  seventy-two  hour  isentropic 
sajnple  was  too  small.  In  each  case  tested,  the  distribution  was 
shown  to  be  reasonably  from  a  normal  population.  In  cases  of  data 
from  isentropic  charts,  there  was  little  doubt  that  the  distributions 
were  from  noriiial  populations.  But  in  cases  from  the  500  millibar 
charts,  the  values  of^  were  close  to  the  critical  values  and, 
therefore,  there  is  some  doubt  that  they  may  be  from  a  normal  popu- 
lation. However,  since  the  values  of  V  did  fall  below  the  critical 
values,  we  will  treat  such  data  as  being  from  normal  populations. 

The  "F"  testis,  pp  152-154J  is  a  statistical  test  for  comparing 
variances  (S  )  of  two  samples  of  data  to  determine  if  the  two   samples 
can  reasonably  be  from  the  same  normal  population.  The  basic  require- 
m.ent  of  the  "F"  test  is  that  the  data  samples  must  each  be  from  normal 
populations,  although  not  necessarily  the  same  norroal  population.  This 
is  why  it  was  necessary  to  show  normality  of  the  data  samples  above. 
To  apply  the  "F"  test,  we  first  postulate  that  there  is  no  significant 
difference  between  samples  tested,  i.e.,  5,^=5^  .  This  is  the  null 
hypothesis.  If  the  value  of  F  computed  falls  below  a  critical  value, 
the  null  hypothesis  is  not  denied  and  there  is  no  statistical  difference 
between  the  two  samples.  If  the  value  of  F  computed  is  higher  than  the 
critical  value,  the  null  hypothesis  is  denied,  and  there  is  a  statisti- 
cally significant  difference  bet\ireen  the  two  samples  and  they  are  not 
reasonably  from  the  same  parent  noniml  population.  It  is  this  last 
condition  we  wish  to  obtain  here.  The  results  of  applying  the  "F'» 
test  are  as  follows: 
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1.  Twenty-four  Hour  V/ind  Direction  Error.  The  varismce  of 
data  from  the  isentropic  charts  5,  -^  /.I.4  .  The  variance 
of  data  fror.i  the  500  millibar  charts  5^'*'=fe.37  .   then 

^  .,;;  a  95So  level  of  belief  is  4.90.  Therefore,  the 
samples  tested  are  significantly  different  and  the  null 
hypothesis  is  denied. 

2.  Forty-eight  Hour  Wind  Direction  Error.  The  variance  of 
date  fror.i  the  isentropic  charts  5,'^==   3.  I  5   .  The 
variance  of  data  from  the  500  millibar  charts  5^\r30- 12.  . 
Then 

f^  at  a  95;o  level  of  belief  is  3.75.  This  is  higlily 
significcint  and  the  null  hypothesis  is  certainly  denied. 

3.  T\-7enty-four  Hour  V/ind  Speed  Error.  The  variance  of  data 
fror.i  the  isentropic  charts  is  155.7.  The  variance  of  data 
fror.  the  500  roillibcir  charts  is  413.1.  Then 

F-    2X5. 

F<„  at  a  95;^  level  of  belief  is  5.19  and,  therefore,  the 
null  hj'pothesis  is  not  denied  and  there  is  no  significant 
difference  bet'.veen  the  saiaples  of  data. 

4.  Forty-eight  Hour  ■■/ind  Speed  Error.  The  variance  of  data 

from  the  isentropic  charts  is  127.3.  The  variance  of  data 

from  the  500  mllibar  charts  is  515.5.  Then 

F=  4.05. 
At  a  90;5  level  of  belief  the  value  of  F\— 4.00  ;  and  at 

this  level  the  difference  between  the  saiaples  becomes 

significant. 
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Significant  results  were  obtained  in  the  cases  of  wind  direction 
error  and  in  one  case  of  \dLnd  speed  error.  In  each  of  the  cases 
where  significance  was  obtained,  an  intuitive  examination  of  the 
data  shows  that  it  is  the  data  from  trajectories  constructed  on 
isentropic  charts  that  are  better. 

It  is  gratifying  to  note  that  the  greatest  iir^provenent  was 
shoMi  for  vrhtd   direction  forecasting  for  a  fortj'"-eight  hour  period. 
In  much  operational  forecasting,  wind  direction  is  usually  deemed 
more  important  than  wind  speed  for  prognostic  work.  The  proper 
placing  of  troughs  and  ridges  and  associated  weather  patterns 
depends  greatly  on  correct  wind  direction  forecasts.  Therefore, 
from  the  results  obtained,  the  CAV  trajectory  method  will  show 
greatest  improvement  over  the  present  system  of  taking  CAV  tra- 
jectories on  constant  pressure  charts  if  used  on  isentropic  charts 
to  forecast  wind  direction,  at  least  up  to  forty-eight  hours. 

A  recent  study  of  CAV  trajectories  along  a  psuedo-600  irdllibar 
chart  was  published  by  Bruch[_?j.  Bruch^s  results  on  wind  direction 
error  are  quite  siirjLlar  to  results  obtained  herein  for  the  500  millibar 
chart.  3ruch  did  not  verify  wind  speeds. 

It  is  realized  that  the  isentropic  chart  leaves  much  to  be  desired 
as  an  operational  ^feather  chart.  At  present  the  chart  is  of  course 
laborious  to  plot  since  no  isentropic  data  are  transrJLtted  on  tele- 
type circuits;  however,  if  such  data  were  transmitted  on  teletype 
circuits,  streai;iline  analj'-sis  vrould  take  little  longer  to  prepare 
than  the  present  500  millibar  chart.  Another  difficulty  is  that  the 
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novenent  of  the  isentropic  surface  is  not  easy  to  predict.  Since 
the  CAV  trajectory  developed  herein  is  a  two   dimensional  projection 
of  a  tiiree  dimensional  path,  the  height  of  the  parcel  in  the  future 
is  in  doubt.  Of  course,  if  the  isentropic  siirface  did  not  move,  the 
future  position  of  a  parcel  would  be  certain.  But,  this  is  not 
generally  the  case.  A  future  study  to  correlate  this  study  with  the 
prognosis  of  isentropic  surfaces  could  be  to  advantage.  Nonetheless, 
the  computation  of  GAV  trajectories  on  an  isentropic  surface  can  be 
of  value  to  prognosis  of  a  nearby  constant  pressure  siu'face.  The 
positions  of  prognostic  troughs  and  ridges  of  the  isentropic  flow 
can  be  placed  on  a  nearby  constant  pressiu'e  surface. 

From  study  of  the  CAV  trajectories  on  the  isentropic  charts  of 
this  paper,  it  qualitatively  appears  that  parcels  ir.ibedded  in  north- 
westerly flow  descend  as  they  move  southeastward,  whereas  parcels 
imbedded  in  southwesterly  flow  rise  as  they  move  northeastward. 

In  the  final  forr.i  of  the  vorticity  equation  there  is  a  diver- 
gence term  that  is  assumed  to  be  zero.  This  assuription  is  not  always 
fulfilled  in  the  atmosphere.  The  effect  of  this  divergence  term  may 
be  important.  For  illustration,  let  the  divergence  ten-'  ^  .^ 
be  constant  in  equation  (9)  and  integrating  xidth  respect  to -fc  ,  we 
obtain 
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Let  "^rrO  at  the  initial  inflection  point,  where  also  P-q  .   'hen 


and 


As  can  be  seen,  positive  values  of  the  divergence  temi  in  a  north- 
westerly flow  will  caiise  a  negative  contribution  to  relative 
vorticity  of  a  parcel  causing  it  to  go  farther  south  than  under  a 
CAV  trajectory.  Positive  values  of  the  term  in  southwesterly  flow 
will  have  a  similar  effect  on  the  relative  vorticity  of  a  parcel 
causing  it  to  not  reach  the  ma>djauni  latitude  indicated  by  a  CAV 
trajectory.  Negative  values  of  the  divergence  term  will  cause 
positive  contributions  to  the  relative  vorticity  of  a  parcel.  In 
northwesterly  flow  this  \d.ll  cause  a  parcel  to  curve  r.:ore  cycloni- 
cally  than  would  be  indicated  by  a  CAV  trajectory  and  to  curve  less 
cyclonically  in  southwesterly  flo\ir.  Fultz  Va,   pp  32-109J  mentions 
that  in  a  lower  level  (10,000  foot  level)  the  average  type  of 
deviation  from  CAV  trajectories  is  such  that  could  be  caused  by 
horizontal  velocity  divergence  in  northwesterly  flow  and  convergence 
in  southwesterly  flow.  Now, 

*  OT     ^X     [c>i.    o^j    ycieix   c>i  c>T 

Thus,  the  divergence  temi  is  composed  of  two  parts,  horizontal 
velocity  divergence  and  a  shear  temi.  Fleagle  Qj  and  Panof  sky  [^TJ 
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have  given  several  values  for  horizontal  velocity  divergence, 
with  a  range  of  about i^  5  x  10"^  per  second.  The  values  of  the 
shear  tern  could  be  computed  at  initial  inflection  points  for 
thirty-four  of  the  CAV  trajectories  on  isentropic  charts  of  this 
study.  In  fourteen  cases  the  ter)i\  was  less  than  10"^  per  second. 
In  eleven  cases  the  value  of  the  shear  term  was  approxiriately  10"^, 
and  in  nine  cases,  approximately  -10"^  per  second.  Usually,  the 
shear  terrr  was  positive  or  negative  if  the  flow  was  up  or  down  the 
isentropic  surface  respectively.  From  this  we  see  that  in  about 
60fo  of  the  cases  studied,  the  shear  terrj  was  of  the  sajne  order  of 
magnitude  as  the  horizontal  velocity  divergence.  In  these  cases 
the  shear  term  either  canceled  or  doubled  the  horizontal  velocity 
divergence  term:.  Now,  if  the  assumption  of  \|!^.^r=-0  were   the 
most  important  assumption  in  the  developraent  of  CAV  trajectories,  we 
might  expect  forecasting  of  future  positions  of  parcels  by  CAV  tra- 
jectories on  isentropic  surfaces  to  give,  in  som.e  cases,  greater 
error  than  from  applying  CAV  trajectories  to  constant  pressure  sur- 
faces. This  was  not  observed  in  the  data  collected  in  this  study. 
This  may  suggest  that  the  assumption  of  no  horizontal  divergence  is 
not  too  restrictive,  and  that  possibly  one  of  the  other  assumptions 
made  by  Rossby  is  more  restrictive,  e.g.,  assumption  of  purely 
horizontal  motion. 

If  we  take  a  large  constant  value  for  horizontal  velocity 
divergence,  say  5  x  10~"  per  second,  take  a  like  value  for  the 
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shear  term,  assume  a  tine  period  of  t\7elve  hours,  and  use  these 
values  in  equation  (11),  we  would  find  that  the  value  of  -f^    at, 
say,  45°  is  reduced  by  33;o  or  to  the  value  of -f  at  latitude  33®. 
But,  this  is  an  extreme  case,  and  the  effects  of  the  divergence 
term  are  usually  m.uch  less.  However,  since  the  assumption  of 
Vj©'^— O  does  lead  to  error,  it  seems  that  some  future  re- 
search could  be  profitably  spent  introducing  quantitative  values 
for  the  divergence  term  and  computation  of  a  trajectory  incorpor- 
ating this  correction  to  the  CAV  trajectory  for  use  in  more 
accurate  forecasting. 
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Figure  2. 
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Figure  3. 
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